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Abstract. We prove that there exists a pair of non-isospectral ID semiclassical Schrodinger 
operators whose spectra agree up to 0{h°°). In particular, all their semiclassical trace in- 
' I variants are the same. Our proof is based on an idea of Fulling-Kuchment and Hadamard's 

variational formula applied to suitable perturbations of the harmonic oscillator. 

(N 

^ 1. Introduction 

C/^ This paper concerns the ID semiclassical Schrodinger operator 

*g Py = -h^-^_ + V(x), ft > 0, 



s 



dx 

where the potential V (which is always assumed to be independent of h) satisfies: 



^ 1/eC°°(M;M), limV^(x) = oo. 

\Q For any h > 0, the spectrum of Py on M is discrete and simple, and we write it as 

a\ 

O spec(Pv') = Ai < A2 < A3 < ^00. 

^ Each Xj depends on h, but we do not include this in the notation. We denote by Uj the 

I corresponding eigenfunctions (which also depend on h), so that 

PyUj = XjUj, Uj E L^(]R). 

^ Our main result, which was conjectured by Colin de Verdiere in jColOSj . is: 

Theorem. There exists a pair of potentials V (x) G C°°(]R) with V (x) > such that the 
operators 

Py.=-h^^^+VHx) (1.1) 

satisfy "spec{Pv+) = spec{Pv-) up to 0{h°°)" (so in particular they have the same semiclas- 
sical trace invariants) and such that the ground state eigenvalues A^ and are different for 
all h > except possibly for a sequence hk 0. 
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Here by "spec(Py+) = spec(Pv'-) up to 0{h°°)" we mean that for every E > and > 
there exists a constant C such that 

sup |A+-A-| <C/i^. (1.2) 

The inverse spectral problem for ID semiclassical Schrodinger operators asks whether spec(Pv') 
determines V uniquely up to a translation x —> x — xq or reflection x — )■ —x. The main tools 
in studying inverse spectral problems are trace invariants such as heat, wave or Schrodinger 
trace invariants. We recommend the surveys |Zel04j by Zelditch, and |DaHell] by Datchev 
and the second author for applications of different kinds of trace formulas in inverse spectral 
results. However the theorem above shows the limitations of semiclassical trace invariants 
meaning that two Schrodinger operators can have the same semiclassical invariants but have 
different spectra. 





Figure 1. Two domains Q and Q' with Tr(cos(tv^A^)) - Tr(cos(tA/A^)) e 
C°^(M), but spec(An) ^ spec(An'). 

An analogous result in the case of bounded plane domains was proved by Fulling and Kuch- 
ment in |FuKu05] where they find two bounded plane domains (called Penrose-Lifshits mush- 
rooms) for which all the wave trace invariants agree but the ground state eigenvalues dis- 
agree. A Penrose-Lifshits domain f2 is a semi-ellipse with asymmetrical bumps A, B and C 
attached to its boundary as in Figure [Tj If one detaches the bump C, reflects it about the 
axis of the ellipse and reattaches it, one gets a non-congruent domain Vt' which has the same 
heat and wave trace invariants as f2, i.e. is indistinguishable from f2 by standard inverse 
spectral techniques. In |FuKu05] . Fulling and Kuchment prove Zelditch's conjecture (see 
[Zel04j ) that these domains are not isospectral. In |Col08j . Colin de Verdiere constructs 
an analogue of the Penrose-Lifshits example for the ID semiclassical Schrodinger operator 
(see Figure [2]). Namely he attaches two small bump functions to the harmonic oscillator 
potential Vq{x) = x"^. By detaching one of these functions, reflecting it about the y-axis and 
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then reattaching it, he is able to construct two non-isomorphic potentials that are isospectral 
modulo h°°. Our goal in this paper is to prove a FuUing-Kuchment theorem for this example 
and verify Colin de Verdier's conjecture that these "isospectral modulo h°°" potentials are 
not isospectral. 




Figure 2. Two potentials with spec(Pv'+) = spec(Py_) up to 0{h°°) but 
spec(Py+) 7^ spec(Py_) 

To our knowledge there are no counterexamples to spectral uniqueness of smooth semiclassi- 
cal Schrodinger operators. There are in fact many recent positive results in this area which 
use trace invariants but make strong assumptions of analyticity or symmetry on the poten- 
tial. In |GuUr07] . the first author and Uribe show that any real analytic potential V on 
M", symmetric with respect to all coordinate axes and with a unique global minimum, is 
determined within the class of all such potentials by the spectrum at the bottom of the well 
of its semiclassical Schrodinger operator. In |Hez09j . the second author and in |CoGu08] . 
Colin de Verdiere and the first author remove the symmetry assumption in dimension 1, 
but keep the assumption of analyticity (see also [Hez09t IGuUrll] for the higher dimensional 
case). In |Col08j . Colin de Verdiere removes both the symmetry and analyticity assumptions 
in the one dimensional case, but adds a genericity assumption and uses all eigenvalues in- 
stead of only the low lying ones. In |GuWa09] . the first author and Wang give a new proof of 
Colin de Verdiere's result with slightly different generic conditions. In |DaHeVe"TT] . Datchev, 
Ventura and the second author show that radially symmetric potentials in M" are spectrally 
determined among all smooth potentials. 

1.1. Remarks. We close the introduction by listing some remarks and related problems: 

• Our method cannot eliminate the possibility of the existence of a sequence hj. 
where the ground state eigenvalues agree. We believe that such a sequence does not 
exist and we actually expect a much stronger statement to hold: 

For all j > 1, there exist cj, Cj > : | A+ - Aj | > Cje'"'/^. 

• We can also ask the same question for semiclassical resonances. Can we find two 
smooth compactly supported potentials and V~ where Pv+ and Py- have differ- 
ent resonances but where the resonances agree up to 0{h°^)l 

• It would be interesting to study the analogous problem in the case of compact Rie- 
mannian manifolds. To be more precise, can one find a pair of Riemannian manifolds 
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{M,g) and {M',g') where spec(Ag) 7^ spec(Ag/) but all the wave trace invariants are 
the same, that is to say Tr e"**v^ - Tre"'*VV ^ C°^(M)? 

We are grateful to Daniel Stroock for bringing to our attention a reference for the Kato- 
Rellich theorem. We are also thankful to Kiril Datchev for his useful comments on the 
earlier version of the paper. The second author would also like to Ramis Movassagh for 
helping us to find a graph of the Weber function using Maple. 



2. Proof of Theorem 



Proof. We will show the existence of the potentials in Theorem [T] by choosing suitable 
perturbations of the harmonic oscillator Vq{x) = x^. Suppose a, (3 E C^(R) such that 
supp(a) C (—3, —2) and supp(/3) C (3,4) and they are not identically zero. We put 

V+{x)=x^ + ta{x) + el3{x), (2.1) 

V~{x) = + ta{x) + el3{-x). 

We denote by Xf{e) the ground state eigenvalue of Py±. We will show that there exist a, j3 
and t, e > small enough such that the pair has the properties asserted in the theorem. 
From the beginning we assume that e and t are small enough that do not have any 
critical points beside x = 0. 



We first show that "spec(Py+) = spec(Pv'_) up to 0{h°°)" in the sense of (1.2). This is 
claimed and proved in |Col08j using the method of "Bohr-Sommerfeld quantization condi- 
tions to all orders" . It can also be proved using quantum Birkhoff normal forms and their 
equivalence with trace invariants as in |GuWa09] . In that paper the first author and Wang 
consider the spectral density measure Uh defined by 

00 

Mf)= Tr/(Py) = ^/(A,) 

i=i 

and they show that it satisfies an asymptotic expansion of the form: 

00 

^/,(/)~(27r/^)-"^afc(/)/^2^ 

fc=0 

where are called semiclassical trace invariants. They show that these invariants are of 
the form: 

where Pkj are universal polynomials. They also give an algorithm to compute a^. Then they 
show that Uh determines the global Birkhoff normal form (whose existence they also prove in 
this setting) of Py up to 0{h°°) and thus determines the eigenvalues up to 0{h°°). However 
for our potentials and V" , the trace invariants in (2.2) are identical and hence their 
are the same up to 0{h°°) which proves the claim. 



A FULLING-KUCHMENT THEOREM FOR THE ID HARMONIC OSCILLATOR 



5 



To show that the ground state eigenvalues A^(e) and (e) are different we use Hadamard's 
variational formula (Lemma 2.3 below) which implies that 

d 



de 



Af(e)= / P{±x)Mx)\'dx, 



(2.3) 



where ui{x) is an normalized eigenfunction of —h"^-^ + + ta{x) with the ground state 
eigenvalue Ai 
enough \ui{x) 
such that: 



Aj'(O) = A;^ (0). Then in Lemma 2.4 we show that for /i = 1 and t small 



^ is not an even function in (—4, —3) U (3,4). Thus by (2.3) we can find a /3 

d 



d 
de 



=0 



de 



=0 



This implies that for e small enough (and t small enough) we have: 

For/i = l: A+(e)^Ar(e). 

Finally by the Kato-Rellich theorem (see |ReeSij Theorem XIL8.) the eigenvalues are 
analytic functions of h for all /i 7^ 0, therefore A^(e) 7^ A]~(e) for all /i > 0, except possibly 
for a sequence /i^ — )■ 0. This finishes the proof of the theorem. □ 



It now remains only to prove Lemmas 2.3 and 2.4 below. Before doing this we makes some 
remarks regarding the proof above. 

Remark 2.1. We note that the Kato-Rellich theorem is very important in our proof. We 
could not follow our argument with a variable h because e would depend on h. We choose 
and fix an e for h = 1 and then use the analyticity in h to argue that, except for a sequence 
hk — )■ 0, for all h > we have Xf{e) 7^ A|f (e). 

Remark 2.2. We also point out that the eigenvalues Xf{e) are not analytic at h = 0. We 
can see this using the theory of quantum Birkhoff normal forms at the bottom of the well of a 
potential which was developed by Sjostrand in |S j92| . Since V^i^x) and Vo(x) = have the 
same Taylor coefficients at a; = 0, the bottom of of their wells, they have the same QBNFs 
at (x,^) = (0,0) and therefore the low lying eigenvalues (in particular the ground states) 
of their Schrodinger operators Py± and Pv^ must have the same asymptotic expansion of 
the form qih + q2h^ + . . . as /i — )■ 0. However the ground state eigenvalue of Vq = x'^ is h 
hence if A^ (e) was analytic at h = then we would have A^ (e) = h. But this is not the case 
if for example we choose a and /3 nonnegative and not identically zero. In fact under this 
assumption we have Xf > h. To see this we recall that 

A+= rnin_ {(-h'f^+x^ + pix))(l),(l))^,. (2.4) 



'L 



2=1 dx"^ 



where p = ta + e(3 and by assumption p > 0. Let ui be an normalized ground state 



eigenfunction of —h'^j^ + + p(x). Then in (2.4) the minimum is attained by ui and 



d'^ f 

The first term is greater than or equal to h with equality only if Ui is a ground state 
eigenfunction of the harmonic oscillator in which case the second term is not zero. So 
A+ > h. 
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We now state Hadamard's variational formula for Schrodinger operators. It is a formula for 
the first variation of the eigenvalues of a perturbed operator in terms of the eigenfunctions 
of the unperturbed operator. Since in this case (unlike for example in the case of bounded 
domains) the argument is simple we give its proof. 

Lemma 2.3. Let V e C"»(M;M) wtth lim\^\^^V{x) = oo. Let (5 e C^{^). Suppose \j{e) 
is the j-th eigenvalue of —h^^ + V{x) + e/3(x) and suppose Uj{x) is an normalized 
eigenfunction of —h'^-^ + V{x) with eigenvalue Aj(0). Then 




Proof. Since in dimension one the eigenvalues are simple, for a given j we can choose a smooth 
one parameter family Uj{x, e) of normalized real eigenfunctions of —h^^ + V{x) + e/3(a;) 
with eigenvalues Aj(e). So by our notation Uj{x) = Uj{x,0). We now write 



d 
de 



A,(e) 



e=0 



d 

de 



e=0 



((-/^'^ + V{x) + eP{x))uj{x, e), u,(x, e)) 



L2 



(/3(x)mj(x),Mj(x))^2 + {{-h 



dx"^ 



Vix))i 



de 



Uj{x,e)),Uj{x)) 



L2 



e=0 



Uj{x,e)) 



e=0 



L2 



Because V is real valued the operator —h"^-^ + V{x) is symmetric and therefore the last two 

^Uj{x,e))^, = 
1. □ 



terms are identical. In fact each one is zero. This follows from [uj{x) 
which in turn follows by applying to the equation {uj{x, e), Uj{x, e)) 



' de I 



'L2 



In the next lemma we put h = 1. 

Lemma 2.4. Let a G C(f (— 3, —2) be nonnegative and not identically zero and let V{x) = 
x^ + ta{x). Suppose ui{x) is a ground state eigenfunction of —-^ + V{x) with the ground 
state eigenvalue Ai. Then there exists t > small enough such that \ui{x)\'^ is not an even 
function on (—4, —3) U (3,4). 



Proof. First of all since away from supp(a) the potential V{x) is real analytic, if |Mi(a;)p 
is even on (—4, —3) U (3,4) then by analytic continuation it is even on (— oo, —3) U (3, oo). 
We also note that (^1(0;))^ being even implies Ui{x) is even. This is because a ground state 
eigenfunction never vanishes. Therefore Ui{x) can not change sign and in particular can not 
be odd anywhere. 

To prove ui{x) is not even on (—00, — 3)U(3, 00) we introduce the parabolic cylinder functions 
below and review some of their basic properties in ^ 



Let W{x) be the unique solution to 

- W"{x) + (x^ - Xi)W{x) = 0, 



(2.5) 
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with 



\imW{x) = and W (-3) = ui{-3) . 

x—^—oo 



The function W{x) is called a Weber function. An argument using (2.4), similar to that 
in Remark |2.2 , shows that for t small enough we have 1 < Ai < 3. In fact we can make 
Ai arbitrary close to 1. By the WKB method we see that |H^(x)| either grows or decays 
exponentially as x — )■ oo, but exponential decay is ruled out by the fact that it would make 
W an eigenfunction of the harmonic oscillator with eigenvalue Ai {1,3,5,...}. Thus 
W{x) is exponentially decaying near — oo with lima;_j._oo W = and exponentially growing 
(in absolute value) near oo. In fact we will see in ^that lim2^_5.oo = —oo, and that W 
has a unique critical point (which is a global maximum) at x = —a, \a\ < \/Xi, and vanishes 
only once for a large positive value of x. In particular, W'{x) is positive for x < —a and 
negative for x > —a. See Figure [3] for a graph of W. We will prove these claims about W in 



Since ui satisfies the same equation as W on {—oo, —3) U (—2, oo), using u{—3) = W{—3) 
we have 

ui{x) = W{x), X < —3, 

and using limj._j.oo Ui{x) = lim^^oo W{—x) = we have 

ui{x) = cW{-x), X > -2 

for a positive constant c. In particular Mi(0) = cW{0). We claim that c > 1 which will show 
that Ml is not even on (— oo, —3) U (3, oo), completing the proof of the lemma. 



The proof that c > 1 uses Sturm-Liouville theory. Let Q{x) 
Xi — x^ — ta{x). Then 

W" + QW = 0, 

and 

u'l + QiUi = 0. 



Ai — and Qi{x) 



We rewrite (2.6) and (2.7) as a first order system using the Priifer substitution: 



W{x) 
W'(x) 



r[x) smt^(xj 
= r(x) cos 9{x) 



Ui{x) = ri{x) sin6'i(x) 
u[{x) = ri(x) cos 6^1 (x) 



(2.6) 
(2.7) 

(2.8) 



where we choose the branches of 6 and 6i such that ^^(— 3) = ^i(— 3) G (0,7r/2) (this is 
possible because W{—3) > and W'{—3) > 0, see Then (see |BirRot] section 5, chapter 
10): 

e'{x) = Q(x) sin^ e{x) + cos2 ^(x), 
6[{x) = Qi{x) sin^ 6i{x) + cos^ Oi{x), 



(2.9) 



and by (2.8) 



(2.10) 



W'{x) = cote{x)W{x), 
u[{x) = cot 6'i(x)ui(x). 

Because Qi{x) < Q{x) and because ^i(— 3) = 6'(— 3), by applying an elementary comparison 
theorem to (2.9) (see |BirRott Theorem 7, Chapter 1]) we get 

Oiix) < e{x) for X > -3. (2.11) 
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< ^ for 



Moreover since W{x) and W'{x) are positive on x < —a (see [|3| we have < 9{x) < | 
X < —a and thus 9i{x) < | for — 3 < s < —a. In fact, because ui never vanishes, 9i{x) is 
never zero and 

< OAx) <- for - 3 < X < -a. 



Using this and (2.11), by applying an elementary comparison theorem to (2.10) and because 
Mi(— 3) = W{—3) we obtain 



Ui{x) > W{x) for 



3 < X < 



-a. 



(2.12) 



This is the key point: (2.12) is an inequality for Ui and W in the region where they solve 
different equations. On the other hand Ui{x) = cW{—x) for x > —2 so 

u\{x) = -cW'i 

-cW'(a) 



-X] 



for X > —2, 



(2.13) 



and in particular u[{—a) = —cW'{a). If 9i{—a) = | then u[ 
W'{a) = 0. But W has only one critical point which is at x : 



= by (2.8) and therefore 
Hence 6'i(— a) < |. Then 



by (2.8) u'{—a) > and thus W'{a) < 0. This shows that a must be positive because if 
a < then the point a is on the left of the critical point —a of W but W is increasing for 
X < —a. Finally W is decreasing on (— a,0) and by (2.13) Ui is increasing on (— a,0). So by 
(|2l2| Mi(0) > W{0) and therefore c> 1. □ 



3. Parabolic cylinder functions (Weber functions) with small frequency 



Let W{x) be a Weber function with small frequency which was defined in (2.5). Here we 



prove the properties of W which are needed in the proof of Lemma 2.4 



That W{x) decays exponentially near — oo and grows exponentially near oo follows from the 
classical WKB method for ODEs (see for example |01] ) . In fact the exact decay and growth 
rates (see |WhiWaj . section 16.5.) are given by 

W{x) ~ C(y2|x|)(^^-^)/2g-xV2 as X ^ -oo, 

W{x)r^ ^.^"^ ^ ry2x)-("^+^)/V'/^ as x^cx), 
1(2(1 - Ai)) 

where C is a positive constant which depends on our normalization W{—3) = ui{—3). Note 
that for Ai > 1 very close to 1 we have r(^(l — Ai)) < which shows that lim^.^00 W{x) = 
—00. 

To prove the other properties of W we first show that W is positive on (—00, 3]. To do this 
we choose rj e C^(3,4) nonnegative and consider the operator 

We represent its ground state eigenvalue by /ii((5) and an associated smooth family of 
normalized eigenf unctions by iIji{x,6). Then by Hadamard's variational formula 

d 



^^-/^i(^) = / 'n{xMi{x,s)\ dx. 
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Figure 3. Weber function with small frequency. 



Similarly 

-Ai(t) = / a{x)\ui{x,t)\'^dx 



dt 

This implies that fii{6) and Xi{t) are increasing functions in 6 and t respectively. Thus 
because /ii(0) = Ai(0) = 1, we can find 6 and t positive so that 1 < fii{S) = Xiit) < 2. 
This is the value t that we choose in Lemma |2.4 Because ipi satisfies the same equation as 
W for X < 3 we have ipi = kW on (— oo, 3]. On the other hand since ipi is a ground state 
eigenfunction it does not vanish on M so the Weber function W has no zeros in (— cxd, 3] and 
must be positive there. 



We can use this fact to show that there is only one critical point which is at the maximum 
X = —a. First of all in the forbidden region x < —^/X we have W{x) > and therefore 
by (2.5) we have W" > there. This implies that W will increase on {—oo, —y/Xi) and in 



particular will never vanish. The maximum x = —a must be in the classical region \x\ < a/Ai 
where W" and W have different signs. There can not be two critical points in the classical 
region because then by Rolle's theorem W" would vanish in between them, and W would 
vanish there as well by the equation, contradicting W > 0. Furthermore, W can not have any 
critical points in the forbidden region x > y/Xi- Assume Xq is such a critical point. Then 
W"{xo) 7^ because if it were zero then by (2.5) W{xo) = which together with would 
W'{xo) = imply that is a trivial solution. If W"{xo) > then by ( [2^ W{xo) > and 
the graph of W should stay concave up for all x > xq which implies that lima;_5.oo W{x) = oo 
which is a contradiction. Finally if W"{xo) < then W^xq) < but because W{3) > we 
should have a zero of W between 3 and xq which shows that xq > 3. Now if zq is the smallest 
zero with Zq > 3 then W'{zq) < 0. After this point Zq, W can only decrease (and is never 
zero) because W" < in this region. This contradicts W'{xo) = 0. 



This concludes the proof of the properties of W which are needed for Lemma 2.4 In the 
course of the proof of this lemma we also prove that a > 0. 
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